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ps) "Abstract 

[We describe an experimental setup specifically designed for measuring the ultrasonic transmission through liquid foams, 
P\J -over a broad range of frequencies (60-600 kHz) . The question of determining the ultrasonic properties of the foam (density, 

[phase velocity and attenuation) from the transmission measurements is addressed. An inversion method is proposed, 

-tested on synthetic data, and applied to a liquid foam at different times during the coarsening. The ultrasonic velocity 
Vw* [and attenuation are found to be very sensitive to the foam bubble sizes, suggesting that a spectroscopy technique could 

-be developed for liquid foams. 
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^p^l. Introduction 
C/2 ■ 

^ [ Acoustic spectroscopy [H Q is a powerful tool to deter- 
mine the particle size distribution in a dispersed media. It 



[consists in measuring the effective attenuation and velocity 
^ -of sound in the medium, over a broad range of frequencies 
•^[(typically 1-100 MHz), and fitting those quantities with 
^!^'the appropriate model in order to obtain the particle size 
[distribution. This technique gives real time measurements 
C^and is not destructive. It does not require dilution of the 
[sample, and recent developments suggest that it could also 
^-H -give information on the local rhcology of the medium Q . 
J> ' Acoustic spectroscopy is routinely used for dispersions 
■of liquid droplets or solid particles, and commercial appa- 



'ratuses are on the market for these kinds of media. Even 



ly-j -gas bubbles can be inspected when their concentration is 
- , • [not too high 0, [HI . But there is no spectroscopy technique 
^~"^ -available for liquid foams, even though those materials 
[represent an important challenge in an increasing number 
y—{ -of industrial areas (food, cosmetic or petroleum industry). 
^ [Absence of ultrasonic spectroscopy can probably be 
• ^ -explained by two difficulties that arise when dealing with 
[the acoustics of liquid foams. First, the models used for 
^ -other media arc not directly applicable __to liquid foams. 



"Specific models have been developed l6|,|7|, but they still 
need to be experimentally validated. Which leads to the 
second problem: acoustic measurements in liquid foams 
are not straightforward. In particular, the measurements 
that exist in the literature 12 | do not cover a 

range of frequencies broad enough to measure a dispersion 
of the acoustical properties, the key feature for extracting 
the size distribution. 

In this article, we propose a new experimental setup for 
measuring the effective velocity and attenuation of ultra- 



sound in liquid foams from 60 to 600 kHz. 



2. Experimental setup 

2.1. Ultrasonic setup 

Our setup is based on the transmission of acoustic pulses 
through the sample to be characterized. As liquid foams 
are mainly composed of air, their acoustic impedance is 
low. Usual ultrasonic techniques that use piezoelectric 
transducers, and water as a coupling medium, face an 
impedance mismatch issue when dealing with such me- 
dia. Air-borne acoustic waves thus seem more appropriate 
in terms of the amount of acoustic energy one can trans- 
mit through the sample. Moreover, the new generation of 
micromachined capacitance air transducers [13[ does not 
suffer from the bandwidth limitation that has made pulse 
techniques impossible for air-borne ultrasound. An ad- 
ditional benefit of working with air as coupling medium 
is to suppress the problem of water infiltration inside the 
sample, an issue particularly important for liquid foams. 

Figure [T] shows a scheme of the ultrasonic setup. A func- 
tion waveform generator, controlled by a computer, gener- 
ated a gaussian pulse, amplified by a power amplifier (up 
to 150 V for the most attenuating samples) and transmit- 
ted to the emitting air transducer (BAT-1, MicroAcous- 
tic). The ultrasonic pulse then propagated through the 
sample before being received by an other transducer, pre- 
amplified (40 dB) and recorded by a digital oscilloscope 
connected to the computer. Two gaussian pulses, cen- 
tered around 150 kHz and 400 kHz, allowed us to explore 
the complex transmission from 60 to 600 kHz. Note that 
we did not use the total bandwidth of the transducers be- 
cause the attenuation was usually too high for frequencies 
higher than 600 kHz. Distance L between the transducers 
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and the cell was 9 cm, so that the cell was in the far-field 
of the transducers (the diameter of the active surface of 
the transducers was 1 cm). 
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Figure 1: Scheme of the ultrasonic setup. A pair of broadband air- 
transducers measures the complex transmission through the sample 
from 60 to 600 kHz. 



2.1.1. Design of the cell 

Liquid foams are not solid, which means one cannot cut 
a well-defined slab of foam and put it in the ultrasonic 
beam. On the other hand, using solid walls to maintain 
the foam pauses the problem of the impedance mismatch 
between the air and the walls. We designed a specific 
cell with very thin walls, depicted in figure [21 Circular 
holes with diameter D - 1 cm were made in two plastic 
plates (18 X 12cm^, 3mm-thick) and covered with very 
thin PET (polyethylene terephtalatc) films. The thickness 
of the films was chosen as small as possible [h ^ 3 jum) in 
order to make them transparent to air-borne ultrasound. 



were clamped onto the walls to rigidity them during the 
filling process (see figure [U . Then, the liquid foam was 
poured onto one of the walls, two spacers with the desired 
thickness d were placed on the edges, and the cell was 
closed with the second wall. The last step consisted in re- 
moving the rigid parts, making the cell ready for ultrasonic 
measurements. Note that the cell was placed horizontally 
to avoid a gradient of liquid fraction perpendicular to the 
ultrasonic beam. As liquid foams were found to be very at- 
tenuating for ultrasound, spacers were thin: 1.1 +0.02 mm 
or even 0.5 ± 0.02 mm for the most attenuating samples. 
However, it was checked that the diameters of the bubbles 
in the foam were smaller than d (see section [2.2.2p . ensur- 
ing a three-dimensional structure for the foam. Note that 
a benefit of the small thickness was that filling the cell did 
not require a large volume of foam (typically 8 mL were 
enough) . 

2.1.2. Data processing 
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Figure 3: Measured signals through air (reference, solid line) and 
through one PET film (sample, dash line), for a gaussian pulse cen- 
tered at 150 kHz. Inset: Fourier transforms of the two signals. 
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Figure 2: The liquid foam is contained in a thin-wall cell whose 
thickness d is well defined. Two rigid plates are used to rigidify 
the walls when the cell is being filled. A pair of electrodes is also 
included in the cell for electrical conductivity measurements, which 
give access to the liquid content of the foam. 

Filling the cell was a critical task because thin films are 
flexible. Two solid plates that fitted exactly in the holes 



As an illustration of the data processing, we show the 
acquisition of the transmission through one PET film. Fig- 
ure [3] shows the two signals recorded by the oscilloscope 
without (Reference) and with (Sample) the film placed be- 
tween the two transducers, for a gaussian pulse centered 
at 150 kHz. Note that even though the film is not perfectly 
transparent to ultrasound, a significant part of the energy 
is transmitted. The transit time is strongly affected: the 
pulses have a quadrature phase relationship. 

The ratio between the Fourier transforms of the two 
temporal signals gives the complex transmission through 
the sample: T - Tsam/Tref-, where Tref and Tsam are 
respectively the complex transmission without and with 
the sample. Figures and HJa show the amplitude and 
the phase of T as functions of the frequency. Note that 
the time limits of a pulse are not always clear. For in- 
stance, in figure [31 one can wonder if the small bump in 
the signal around 305 /is is part of the pulse, or a spurious 
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Figure 4; Transmission through one film. Amplitude (a) and phase 
(b) of the transmission coefficient are obtained by calculating the 
Fourier transforms of the signals shown in figure [S] Dashed lines 
correspond to the theoretical transmission through a film with thick- 
ness /j = 3.1 yum. 



signal. To solve this issue, we systematically made two 
different truncations: a short one that captured only the 
main pulse, and a broader one that included other small 
features. For the pulses of figurc[31 the first truncation was 
270 to 300 ;us. and the second one from 250 to 320 //s. The 
differences between the two truncations appear as errobars 
in the amplitude and phase of the transmission (very small 
in figures andlUD because the signal over noise ratio was 
good). 

2.2. Foam production and characterization 

As the acoustic properties of liquid foams are known 
to strongly depend on their structure, we included other 
characterization techniques to the setup for measuring pa- 
rameters such as the liquid fraction <D of the foam and 
the bubble size distribution. Liquid foams have been the 
subject of an extensive literature, and many experimen- 
tal techniques for characterizing them have been devel- 
oped The technique used for the production of the 
liquid foam is known to have a great influence on the struc- 



ture of the final product. We used a very simple technique 
known as the "double-syringe" technique. In one syringe, 
one takes a volume V( of liquid and completes with gas 
to a total volume V . Then the second syringe, empty, is 
connected to the first one and the mixture is pushed from 
one syringe to the other several times. With this simple 
process, the gas is sheared and divided into small bubbles. 
The size of the bubbles depend, among other parameters, 
on the size of the aperture that connects the two syringes. 
In our case, the syringes were connected with a standard 
female-female Luer Lock connector; the length and diame- 
ter of the cylindrical shearing zone were 16 mm and 4 mm, 
respectively. Note that the syringes should be rubber free, 
because rubber can limit the foamability of the solution. 

For the foams presented in this article, the liquid phase 
was distilled water with lOg/L of SDS (sodium dodecyl 
sulfate) and 0.5 g/L of xanthane. We used SDS as a sur- 
factant because it gives a good foamability to the liquid, 
and xanthane because it increases the viscosity of the liq- 
uid phase, thus limiting the drainage of the foam [l5| . The 
gaseous phase was air with traces of C(,F\a (perfluorohex- 
ane) to decrease its solubility, which slows down the coars- 
ening of the foam [Tij . 



2.2.1. Liquid fraction measurement 
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Figure 5: Electrical scheme of the conductivity setup. Measuring 
tension U when the cell is filled with liquid and when it is filled with 
foam gives the relative conductivity of the foam, which is related to 
its liquid fraction. 

The advantage of the "double-syringe" technique is the 
relative control it offers on the liquid fraction of the foam, 
as one chooses the amount of liquid and gas taken to make 
the foam (<[) = Vc/V). However, with standard small sy- 
ringes and without specific precautions, the uncertainties 
on the volumes are such that the actual liquid fraction 
can noticeably differ from the expected one. We measured 
the liquid fraction by weighting the foam and found that, 
typically, when we aimed at a liquid fraction of 10%, the 
actual O was between 7 and 13%. In order to have a 
precise and in-situ determination of the liquid fraction we 
set up an electrical conductivity measurement [lil [l7j . As 
shown in figure [21 two electrodes (5 mm in diameter) were 
included in the walls of the cell. They were placed as close 
as possible from the edge of the film (3 mm in practice), 
so that they probed the same part of the sample as the 
one probed by the ultrasonic beam. Figure [S] presents the 
scheme of the electrical setup. The generator is set on an 
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alternative tension of 1 kHz, with a voltage of £ = 1 V and 
the resistance of the circuit is R = 1.2 kQ. The relative 
conductivity of the foam cr,. is given by 

Veil +RiJR) 



Uf_ E_ 

ll ^ E-U{\ +Rin/R) 



(1) 



where Rin is the internal resistance of the generator, U the 
tension measured on the resistor when the foam is in the 
cell, and Uc the tension measured when the cell is filled 
with liquid. From the relative conductivity^ liquid fraction 
(S) is obtained thanks to the empirical law 

3cr,(l + 11 cr,) 
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<5 = 



1 +25cr,+ 10cr2' 



(2) 



For the liquid foam we present in this article, we aimed at a 
liquid fraction of 5 % and obtained O - 6 + 0.5 % according 
to the conductivity and weight measurements. 

2.2.2. Bubble size measurement 
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Figure 6: Image analysis of bubble rafts (see inset) gives the his- 
togram of the bubble radii distribution. Here the distribution is fitted 
with a log-normal law (solid line); the median radius is IS/im and 
the polydispersity is 0.45. 

Measuring the bubble size in a three-dimensional foam 
is very challenging. We adopted a quite intrusive method 
that consists in pouring a small volume of the foam onto 
the surface of the same liquid used for making the foam. 
One can then spread the bubbles in order to obtain a two- 
dimensional bubble raft, easy to image (see inset of fig- 
ure ini) • If it is assumed that the foam that remains in the 
syringe has the same time evolution than the one that has 
been put into the cell, this technique can even allow one to 
follow the evolution of the size distribution with time. Fig- 
ure [6] shows an histogram of the bubble radii, obtained for 
our example foam at f = 0. It can be fitted by a log-normal 
law: 

n,,, I (ln(r/ro))^ \ 

n{r)^ ;=SXP ' (3) 

reV27r \ 2e- / 

where n,ot is the total number of bubbles, ro the median 
radius of the bubble, and e the log-normal standard devia- 
tion. Parameter e is an indication of the polydispersity of 



the foam: e — 0.5 means that there is a significant number 
of bubbles whose radius is half or twice the median radius. 



3. Transmission measurements 

3.1. Ultrasonic characterization of the cell 

Being able to measure the acoustic properties of the 
foam contained in the cell requires a good knowledge of 
the cell itself. The first step is to measure the transmis- 
sion through one film (see figure |4]). As the impedance 
of PET is much larger than the one of air (see table [T]), 
and the thickness of the film is small compared to the 
wavelength in air {A ^ 0.5 mm at 600 kHz), the acoustic 
transmission through the film is well approximated by (see 
appendix [Appendix A ) : 



Tm-n 



-^\(kp~ka)h 



1 



(4) 



where / is the frequency, and k, Z, p the wave numbers, 
acoustic impedances, and mass densities, respectively. In- 
dices a and p refer to air and PET, respectively. As shown 
in figure SI equation (0]) fits the experimental data very 
well for h — 3.1 yum. 



Medium 


index 


phase velocity 


density 


air 


a 


340 m/s 


1.2kg/m^ 


PET 


P 


2540 m/s 


1400 kg/m^ 



Table 1: Physical parameters for air and PET. 

The amplitude of transmission through the film de- 
creases as frequency increases: at 600 kHz, only 6 % of 
the amplitude is transmitted through the film. This is a 
limitation of our setup for high-frequency measurements. 
Note that the phase of the transmission also depends on 
frequency (see fig. IUd) and tends toward 7r/2 at high fre- 
quencies. It means that even though the film is very thin, 
Tmm is not purely real, which is an important feature for 
the inversion procedure (sec section 2]). 

To make sure the roles of the two films were well un- 
derstood, we also measured the transmission through the 
empty cell. Figure [7] shows the amplitude and the phase of 
the transmission through this empty cell with a,d - 1.1 mm 
spacer. Multiple reflections occur and give rise to very 
clear Fabry-Perot resonances. The experimental data 
are successfully compared with the theoretical prediction 
Tcc\\{ka,Pa) (see appendix |Appendix A ), which means that 
the cell is well characterized and ready for being filled with 
foam. Note that the first peak of transmission through 
the empty cell is particular. Indeed, the Fabry-Pcrot cri- 
terion for the resonance to occur can be simply written 
as kad — nn with integer « > 1, which gives frequencies of 
maxima for f„-nx 155 kHz. It does not predict the first 
maximum which occurs at 37 kHz. This first maximum is 
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actually a signature of the non purely real reflection coef- 
ficient of the films, which brings an additional phase shift 
and changes the condition of resonance. 
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Figure 7: Transmission through the empty cell. Amplitude (a) and 
phase (b) of the complex transmission though the empty cell are 
shown. The circles are data points, the continuous lines the the- 
oretical prediction, and the dashed lines the transmission predicted 
through two films with no multiple reflections. 



3.2. Transmission through liquid foams 

When the cell is filled with a liquid foam, the trans- 
mission does not show any Fabry-Pcrot resonances (see 
figure IS]) , which is not a surprise given the high attenua- 
tion of ultrasound in liquid foams. We show results for a 
foam at two different ages, with a constant liquid fraction 
<J> = 6 %. Note that the attenuation was so high that a thin 
cell was used {d — 0.5 mm). Even with this precaution, the 
signal could not be measured above 180 kHz bX t - Omin. 

4. Data analysis 

4.1. Principle 

The analysis of the experimental data is an inversion 
problem: from the measurement of the complex trans- 
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Figure 8: Transmission through liquid foams. Amplitude (a) and 
phase (b) of the transmission through the cell filled with the same 
liquid foam at two different ages are shown. For comparison, the 
transmission predicted through two films with no multiple reflections 
is reported with dashed lines. 



mission through the cell, can we determine the velocity 
c and the attenuation a of ultrasound in the medium? 
This is equivalent to determining the complex wavenum- 
ber k — a>/c + ia/2. As recalled in the appendix, the to- 
tal transmission through the cell depends on the density 
and wavenumber in all the media involved, and on the 
thicknesses of the different layers. As the cell is well char- 
acterized, only 2 quantities are unknown: p and k, i.e. 
the density and the wavenumber for the medium that fills 
the cell. If we note T^^p the experimental transmission 
measured through the cell, the question is then to know 
whether we can find p and k such that Tccii{k,p) — T^xp- 

A first simple approach is to consider that, as there is 
no multiple reflection in the liquid foam, Tceii(k,p) reduces 
*° -^fiim^'**^ because the acoustic wave just goes across 
the two films and the foam. This leads to 



id 



Tfii, 



(5) 
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However, the transmission through the film is modified 
by the presence of the foam. The actual transmission 
through the cell is given by (see appendix [Appendix A[ ) 



Tcoi\ik,p) 



4Zapa)/k 



(6) 



which is much more challenging to inverse. First, both p 
and k are involved in the equation, which means that one 
of them needs to be determined by another measurement. 
Let us assume that we know p. Even with this simplifica- 
tion, k cannot be extracted analytically from equation 
Wc used an iterative process: the first guess ^'"^ is given 
by eq. ([S]) and the iteration proceeds with 



^[«+i] ^ ^ ln( ^"""P 



(7) 



Usually, 10 iterations were enough to obtain a good con- 
vergence. 

^.2. Test on synthetic data 

We calculated Tcci\{k,p) for different k and applied the 
inversion method to the synthetic data obtained. Conver- 
gence of the method was found to be excellent, providing 
that multiple reflections in the cell could be neglected (i.e., 
ad should be large enough). 

Figure [3] shows a typical result for a test medium with 
p = 50 kg/m\ c = (30 + 500/) m/s and o- = (2 + 8/) mnr' , 
where frequency / is in MHz- As a first step, we calculated 
the complex transmission through a fictive cell filled with 
the test medium (sec appendix [Appendix A[ ) . Amplitude 
and phase of the complex transmission through the cell are 
shown in figure (9] A second step consisted in applying the 
inversion method to those synthetic data, which gave back 
the actual velocity and attenuation (figure ITU)) . 

Note that, as for any phase measurements, the phase of 
the transmission is measured modulo 2;r, which can affect 
the value found for the attenuation and the velocity. As 
an example, figure [TOj shows how the inversion is affected 
when a global 2n phase shift is applied to the phase of the 
transmission. The attenuation is only slightly changed, 
but the velocity is significantly lowered. Determining the 
correct 2n shift for the phase has always been a problem 
for ultrasonic velocity measurements [l8|. In many cases, 
the ambiguity is resolved because one knows that the ve- 
locity should be within a reasonable limit, or because in- 
correct 2n shifts lead to an unexpected dispersion. But 
for media such as liquid foams, many options give realis- 
tic results: both results in figure ITOb look plausible, for 
example. The solution to the problem is found by noting 
that the phase of the transmission should tend to zero for 
zero frequency [l^l ■ Then, providing that broadband mea- 
surements are available, one can unambiguously determine 
the correct phase (see figure [Sb for example). However, it 
requires the assumption that the phase has a smooth be- 
havior for frequencies below the lowest experimental fre- 
quency. 
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Figure 9: Synthetic data are generated by calculating the amplitude 
(a) and phase (b) of T^^ii(k,p) for known values of k and p. The 
phase is known modulo 2n. As an example, the calculated phase 
(rectangles) can be affected by a +2n shift (circles), resulting in a 
different result for the inversion (see text). 



4-3. Experimental results for liquid foams 

As an illustration, we show that our technique allows us 
to follow the time evolution of the acoustical properties of 
a foam. 

Applying the inversion method to the real data of fig- 
ure p requires first to determine p. We took Wood's 
law [l9[ (also known as the mixture law), which assumes 
that the effective mass density is the same as the actual 
mass density: p — Op^ + (1 — '^)Pg, where p[ and pg are the 
mass densities of the liquid and gaseous phases, respec- 
tively. For the liquid foam presented in this article, the 
inversion was thus performed considering that the mass 
density was p - 60kg/ni-'. Note that other models exist 
and predict that the effective density can be complex and 
frequency dependent [20| . 

As mentioned earlier, the phase is measured modulo 
In. As can be seen in figure [SJ), the measured phase at 
t — min needs a +2n shift in order to follow our pre- 



6 




0.2 0.4 
Frequency (MHz) 

(a) 



0.6 



350 




0.2 0.4 0.6 

Frequency (MHz) 

(b) 



Figure 10: When the inversion method is applied to the synthetic 
data of figure \^ one can find back the attenuation (a) and phase 
velocity (b) in the test medium. Different results are obtained for a 
zero phase shift (rectangles) or a 2n phase shift (circles). The true 
attenuation and velocity are shown in solid lines. 



Errorbars were estimated from the three sources of uncer- 
tainty: hquid fraction (+0.5 %), spacer thickness (±20 /zm), 
and accuracy of the transmission measurement (see error- 
bars in Fig.|S]for example). Even though the hquid content 
remains the same, the hquid foam shows very different ul- 
trasonic behaviours as it ages. When the bubbles are small, 
attenuation is large and very dispersive, whereas the veloc- 
ity is almost constant at about 40 m/s. Interestingly, when 
the bubbles grow, the attenuation becomes lower and the 
velocity is found to be dispersive, reaching a plateau at 
about 200 m/s for t - 120 min. This dispersive behavior 
has never been observed in liquid foams. It suggests that 
the size of the bubbles is an important parameter for the 
acoustic properties of liquid foams. 
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Figure 11: Phase velocity calculated from the experimental transmis- 
sion of figure O using the inversion method for different "2jt " shift: 
On (triangles), +2n (circles) and +4n (diamonds). 



5. Conclusion 



scription of a "zero phase at zero frequency". Figure 1111 
brings a confirmation that this is the correct phase cor- 
rection. Indeed, the inversion was applied for 0, +27r and 
+4;r shifts. The Ott shift can be ruled out because it gives 
a very dispersive behavior. The +47t shift gives a plausi- 
ble result, but slightly more dispersive than the +27r shift. 
Dispersion is not the only criterion: the value of the phase 
velocity found with the +2;r shift is in a close agreement 
with Wood's prediction: cwood = V100/f/<l)/(l - <l>) for a 
gas at atmospheric pressure and a liquid with density of 
water, where k — I for isothermal transformations (solid 
line in Fig. [TT|) and k = lA for adiabatic transformations 
(dashed line in Fig. [TT|) . Other experiments at lower fre- 
quencies have found a velocity in_liquid foams close to the 
isothermal Wood's velocity 



Figure [12] shows the attenuation and phase velocity ob- 
tained for the liquid foam from to 120 min, i.e. as the 
bubbles were growing from 13yum (see Fig. ^ to 65 /jm. 



We have described a new experimental setup for mea- 
suring the ultrasonic transmission through liquid foams 
over a large range of frequencies (60-600 kHz). The ul- 
trasonic setup was combined with conductivity measure- 
ments and optical observations that allowed us to charac- 
terize the liquid foam (liquid fraction and bubble size dis- 
tribution). We have shown that, providing that the mass 
density of the foam was known, the velocity and atten- 
uation of ultrasound in the foam could be deduced from 
transmission measurements. As an illustration of the tech- 
nique, the acoustic properties of a foam at different ages 
were presented. The ultrasonic attenuation was found to 
be strongly dispersive and affected by coarsening. For the 
foam with small bubbles (at t — Omin with median radius 
of 13jum) the velocity was well predicted by Wood's ap- 
proximation. On the other hand, for larger bubbles, the 
velocity was very dispersive. Thus, there is a clear sig- 
nature of the bubble sizes in the acoustic measurements, 
which is a good indication that ultrasonic spectroscopy can 
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Figure 12: Attenuation and phase velocity calculated from the exper- 
imental transmission of figure\^ using the inversion method. Atten- 
uation and phase velocity for two intermediate ages of the foam are 
also reported. 



be envisaged in liquid foams. Note that an important fea- 
ture of our experimental setup is real-time measurements 
(one acquisition is taken in less than 30 s). This is a ma- 
jor advantage for inspecting media such as liquid foams, 
whose time evolution is usually significant, or even crucial 
for some applications. 

Further studies are needed for determining the exact 
link between the structure of the foam and its ultrasonic 
response. In particular, a systematic comparison of the 
experimental findings with the theoretical models should 
be undertaken. Hopefully, our setup might help for a bet- 
ter understanding of how ultrasound propagates in a liquid 
foam, an heterogeneous medium whose acoustic properties 
are still not clearly known. 
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Appendix A. Acoustic transmission through a 
three-layer system 

When an acoustic wave propagates through a layered 
media, such as the cell considered in this paper, many 
reflected waves are involved. For 5 media (sec figure [XTS]) 
the total transmission is given by [U 



^5- 



+Z2 

'z'" +z^ 



Z^" + Z3 
Z^^ + Z4 
2Zi 



Z1+Z2' 



(A.1) 



where di is the thickness of layer /, Z,- - pico/k,, ki stands 
for the impedance and the wave vector in medium ; respec- 
tively, and Z™ is the input impedance for layer /, which is 
given by 



Zj^j - iZ,tan(/t,^f,) ^ 



for ; > 1 . 



(A.2a) 
(A.2b) 



Z,-iZ'^jtan(M,) 

With Eqs. (|A.1|) and (|A.2I) . one can predict the transmis- 
sion through the cell filled with a medium whose density 
is p and wavenumber k: Tceii(p,k) — Ti^^/explikaid + 2h)], 
where Ti^s is calculated with ki — ks — kg, p\ — ps — Pa, 
kp, pi — Pa - Pp, d2 - di, - h, d^ - d and k^ - k, 
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Figure A. 13; The transmission through 5 successive layers is consid- 
ered. Note that the layers are numbered backward compared to the 
direction of propagation. 

Several simplifications can occur. A first interesting sit- 
uation is when layer 2 is thin {k2d2 <K 1). Then the trans- 
mission from layer 3 to layer 1 simplifies into 



2Z, 



1 



z, +z, 



., , / z,-+z,z, \ 



(A3) 



If the acoustic impedance of the wall is large compared to 
the ones of the media it separates (Z2 » Z3 and Z2 » Zi), 
another simplification occurs: 



^3^1 = 



2Zi 



1 



Z,+Z3 



(A.4) 
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which gives eq. (|4]) if Zi = Z3. 

Another interesting case is when the attenuation in 
medium 3 is such that muhiple reflections can be neglected 
(note that it also works if c/3 is large enough for multiple 
echos to be separated in the time domain), which leads to 
Zj" = Z3. Thus, for a symmetrical configuration (1 =5, 
2-4) with thin walls {k2d2 <^ 1) of high impedance 
(Z2 » Z3 and Z2 » Zi), if multiple reflection can be ne- 
glected in medium 3, eq. (|A.1|) reduces to 

T,^, = X S!f (A,5) 
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